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ABSTRACT 


We  report  on  the  progress  toward  the  development  of  attenuation  models  for  short-period  (12-22  sec)  Rayleigh 
waves  in  Asia  and  surrounding  regions.  These  models  are  defined  by  maps  of  attenuation  coefficients  across  the 
region  of  study  in  the  specified  period  band.  These  models  are  designed  to  calibrate  the  regional  surface-wave 
magnitude  scale  and  to  extend  the  teleseismic  surface-wave  magnitude — body-wave  magnitude  ( Ms-mb ) 
discriminant  to  regional  distances. 

In  order  to  obtain  accurate  attenuation  estimates,  we  must  first  measure  surface-wave  amplitudes  reliably.  Taking 
advantage  of  certain  characteristics  of  Rayleigh  waves,  such  as  the  dispersion  and  the  elliptical  particle  motion,  we 
employed  a  suite  of  techniques  in  making  accurate  fundamental-mode  Rayleigh-wave  amplitude  measurements.  We 
first  analyze  the  dispersion  of  the  surface  wavetrain  using  a  spectrogram.  Based  on  the  characteristics  of  the  data 
dispersion,  we  design  a  phase-matched  filter  by  using  either  a  manually  picked  dispersion  curve,  or  a  model- 
predicted  dispersion  curve,  or  the  dispersion  of  the  data,  and  apply  the  filter  to  the  seismogram.  Intelligent  filtering 
of  the  seismogram  and  windowing  of  the  resulting  cross-correlation  based  on  the  spectrogram  analysis  and  the 
comparison  between  the  spectrum  of  phase-match  filtered  data  and  raw-data  and  source  spectra  effectively  reduces 
amplitude  contaminations  from  surface-wave  higher  modes,  multipathing,  body- wave  energy  and  other  noise 
sources,  and  results  in  reliable  amplitude  measurements  in  many  cases.  We  implemented  these  measuring  techniques 
in  a  graphic-user-interface  tool  called  Surface  Wave  Amplitude  Measurement  Tool  (SWAMTOOL).  Using  the  tool, 
we  collected  and  processed  waveform  data  for  200  earthquakes  occurring  throughout  2003-2006  inside  and  around 
Eurasia.  The  records  from  135  broadband  permanent  and  temporary  stations  were  used. 

After  obtaining  surface-wave  amplitude  measurements,  we  analyzed  the  attenuation  behavior  of  the  amplitudes 
using  source-  and  receiver-specific  terms  calculated  from  the  three-dimensional  (3D)  velocity  model  CUB2  of  the 
region.  Based  on  the  results,  we  removed  amplitudes  that  yielded  negative  average  attenuation  coefficients,  and 
included  an  additional  parameter  in  the  inversion  to  account  for  the  possible  bias  in  Harvard  Centroid  Moment 
Tensor  (CMT)  (Dziewonski  et  al.,  1981)  scalar  moments.  We  used  the  tomographic  inversion  to  obtain 
surface-wave  attenuation-coefficient  maps  from  12  to  22  seconds  for  Asia  and  surrounding  regions.  The  inverted 
attenuation  maps  are  consistent  with  geological  features  of  Asia.  We  observe  low  attenuations  in  stable  regions  such 
as  eastern  Europe,  the  Siberian  platforms,  the  Indian  Shield,  the  Arabian  platform,  the  Yangtze  Craton,  and  others. 
High  attenuation  is  observed  in  tectonically  active  regions  such  as  the  Himalayas,  the  Tian  Shan,  Pamir  and  Zagros 
mountains.  Finally,  we  conducted  the  calibration  of  a  new  Ms  formula  (Russell,  2006)  for  the  same  region  using  3D 
group- velocity  models. 
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OBJECTIVES 


The  objectives  of  this  study  are  1)  to  develop  short-period  (12— 22s),  two-dimensional  (2D)  Rayleigh- wave 
attenuation  maps  for  Asia  and  surrounding  regions  along  with  associated  uncertainty  statistics  through  a 
tomographic  approach,  and  2)  to  calibrate  Russell’s  (2006)  Ms  formula  with  these  maps  for  the  same  region. 

Knowledge  of  the  seismic-energy  loss  during  the  propagation  of  surface  waves  from  the  source  to  receivers  is 
essential  for  the  accurate  estimation  of  the  surface-wave  magnitude  Ms  and  the  seismic  moment  of  the  source.  This 
is  especially  important  for  nuclear-explosion  monitoring,  in  which  Ms  is  used  in  the  most  robust  seismic 
discriminant,  the  Ms-mb  discriminant.  In  order  to  apply  this  discriminant  to  regional-distance  monitoring,  a  modified 
Ms  formula  using  shorter-period  (<  20  s)  surface-wave  amplitudes  is  required  (e.g.,  Marshall  and  Basham,  1970; 
Russell,  2006).  At  regional  distances,  seismic-wave  propagation  is  strongly  influenced  by  the  lateral  heterogeneity 
of  the  crust  and  upper-mantle  material  properties.  Short-period,  2D  surface-wave  attenuation  maps  developed  from 
observed  amplitude  data  used  to  correct  for  the  propagation  effects  in  calculating  Ms  hold  the  potential  to  reduce 
station-magnitude  scatter  and  network-magnitude  bias.  In  this  paper,  we  describe  the  surface- wave  amplitude 
collection  and  measurement,  the  development  of  the  attenuation  models,  and  the  Russell  (2006)  Ms  calibration 
results. 

RESEARCH  ACCOMPLISHED 


Surface- Wave  Data  Collection  and  Amplitude  Measurement 

Surface-wave  amplitudes  could  be  contaminated  by  a  variety  of  sources  such  as  multipathing,  focusing  and 
defocusing,  body  wave,  higher-mode  surface  wave,  and  other  noise  sources.  In  order  to  obtain  accurate  attenuation 
estimates,  we  must  measure  surface-wave  amplitudes  reliably  by  reducing  the  contamination  as  much  as  possible. 
Taking  advantage  of  certain  characteristics  of  Rayleigh  waves,  such  as  the  dispersion  and  the  elliptical  particle 
motion,  we  employed  a  suite  of  techniques  in  making  accurate  fimdamental-mode  Rayleigh-wave  amplitude 
measurements.  We  first  analyze  the  dispersion  of  the  surface  wavetrain  using  a  spectrogram.  Based  on  the 
characteristics  of  the  data  dispersion,  we  design  a  phase-matched  filter  by  using  either  a  manually  picked  dispersion 
curve,  or  a  model-predicted  dispersion  curve,  or  the  dispersion  of  the  data,  and  apply  the  filter  to  the  seismogram. 
Intelligent  filtering  of  the  seismogram  and  windowing  of  the  resulting  cross-correlation  based  on  the  spectrogram 
analysis  and  the  comparison  between  spectrum  of  phase-match  filtered  data  and  raw-data  and  source  spectra 
effectively  reduces  amplitude  contaminations  from  surface-wave  higher  modes,  multipathing,  body-wave  energy 
and  other  noise  sources,  and  results  in  reliable  amplitude  measurements  in  many  cases. 

We  implemented  these  measuring  techniques  in  a  graphic-user- interface  tool  called  Surface  Wave  Amplitude 
Measurement  Tool  (SWAMTOOL).  Figure  1  is  a  computer-screen  snapshot  of  SWAMTOOL.  The  tool  consists  of 
four  panels  showing  the  seismograms  (upper  right),  a  map  (lower  right),  the  data  spectrogram  (upper  left),  and  the 
spectra  of  the  source,  the  data  and  the  noise  (lower  left).  The  surface-wave  segment  is  first  isolated  with  the 
guidance  of  nominal  group-velocity  marks,  shown  as  red  vertical  dashed  lines  in  the  top  figure  of  the  seismogram 
panel,  and  the  dispersion  characteristics  predicted  by  2D  group- velocity  models  (Ritzwoller  and  Levshin,  1998; 
Levshin  and  Ritzwoller,  2003;  Stevens  et  al.,  2001;  Levshin  et  al.,  2003),  depicted  as  open  red  circles  in  the 
spectrogram  panel.  A  phase-matched  filter  is  then  constructed  with  a  dispersion  curve  (red  line  in  the  spectrogram 
panel)  determined  by  analyzing  the  surface-wave  spectrogram  and  the  dispersion  predicted  by  the  group-velocity 
models.  Depending  on  the  characteristics  of  the  data  dispersion,  we  use  either  a  manually  determined  dispersion 
curve,  or  the  dispersion  curve  of  the  data,  or  the  dispersion  curve  predicted  by  the  group-velocity  model  to  construct 
the  phase-matched  filter.  Surface-wave  spectral  amplitudes  (red  line  in  the  spectrum  panel)  are  measured  after  the 
seismogram  is  processed  with  phase-match  filtering  and  windowing.  We  use  theoretical  source  spectrum  (green  line 
in  the  spectrum  panel)  as  a  reference  for  windowing  the  cross-correlation  resulting  from  the  phase-match  filtering 
(middle  figure  in  the  seismogram  panel),  taking  into  account  the  source-depth  uncertainty  (black  lines  bracketing  the 
green  region  in  the  spectrum  panel.)  We  also  calculate  the  backazimuth  of  the  incoming  Rayleigh  wave  and  compare 
it  with  the  great-circle  backazimuth.  The  result  is  displayed  at  the  lower  left  of  the  map  panel.  Other  information 
displayed  in  the  tool  includes  the  average  Q  between  the  source  and  the  station  at  10,  15,  and  20  s  (in  the  spectrum 
panel),  noise  spectrum  (blue  line  in  the  spectrum  panel),  and  source  and  path  information  in  the  map  panel.  Finally, 
the  usable  frequency  band  of  the  measured  spectrum  is  determined  by  the  yellow  region  in  the  spectrum  panel. 
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Figure  1.  A  computer-screen  snapshot  of  SWAMTOOL.  It  shows  an  example  where  the  medium  between  the 
source  and  the  receiver  is  relatively  simple  with  weak  lateral  heterogeneities.  The  seismogram 
possesses  well-behaved  dispersion  characteristics,  which  are  also  predicted  by  the  group-velocity 
model.  The  spectrum  of  phase-match  filtered  data  does  not  differ  significantly  from  the  spectrum  of 
the  raw  data. 

Figure  1  shows  an  example  where  the  medium  between  the  source  and  the  receiver  is  relatively  homogeneous. 
Figure  2  gives  an  example  illustrating  how  multipathing  is  treated.  An  example  showing  the  reduction  of  noise  from 
body  waves  and  other  sources  using  the  tool  is  given  in  Figure  3. 

We  collected  and  processed  waveform  data  for  200  earthquakes  occurring  throughout  2003-2006  in  and  around 
Eurasia.  The  magnitudes  of  these  events  range  from  5  to  6.  Source  depths  are  less  than  70  km.  Data  from  135 
broadband  permanent  and  temporary  stations  were  used.  Using  SWAMTOOL,  we  made  both  two-station  amplitude 


Figure  2.  This  example  illustrates  how  the  effects  of  multipathing  and  focusing  are  mitigated  by  using 
SWAMTOOL.  In  this  example,  the  surface  wave  from  the  source  traverses  the  Tarim  Basin, 
resulting  in  surface-wave  packets  traveling  along  different  paths  and  arriving  at  the  receiver  at 
slightly  different  times.  The  cross-correlation  from  phase-match  filtering  shows  two  peaks 
corresponding  to  the  two  surface-wave  packets,  although  they  are  difficult  to  distinguish  in  the 
spectrogram  due  to  the  similarity  of  their  travel  times.  Windowing  of  the  cross-correlation  removes 
one  of  the  surface-wave  packets  and  reduces  the  effects  of  multipathing  and  focusing. 
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Figure  3.  This  figure  is  an  example  in  which  the  surface-wave  signal  is  contaminated  by  body  waves  and 

signals  from  other  sources.  The  spectrogram  shows  strong  signals  from  other  sources,  in  addition  to 
the  fundamental-mode  surface  wave.  Because  signals  from  sources  other  than  the  surface  wave  have 
different  dispersion  characteristics,  the  application  of  phase-matched  filter  and  windowing  can 
effectively  remove  these  signals  partly  because  they  appear  as  signal  packets  in  the  cross-correlation 
at  different  lag-times  from  that  of  the  surface  wave. 

ratio  and  single-station  spectral-amplitude  measurements,  but  only  single-station  amplitudes  are  used  in  later 
tomographic  inversions.  Figure  4  shows  the  path  coverage  of  single-station  amplitude  measurements  at  different 
periods. 

Amplitude  Selection  Based  on  Average  Attenuation  Analysis 

Before  using  measured  amplitudes  in  the  tomographic  inversion,  we  estimated  the  average  attenuation  between 
sources  and  receivers  from  the  measurements.  Measurements  that  yielded  negative  attenuation  coefficients  were 
rejected.  We  used  the  theoretical  surface-wave  amplitude  formula  for  a  heterogeneous  Earth  (Woodhouse,  1974; 
Levshin,  1985;  Levshin  et  al.,  1989)  to  estimate  the  average  attenuation.  For  a  laterally  heterogeneous  medium,  the 
spectral  amplitude  of  a  surface  wave  can  be  expressed  in  an  asymptotic  form  as 

A(co)  =  S(co ,  h,  (p)P(a>)B(co),  (1) 

where  co  is  circular  frequency;  h  is  source  depth  and  cp  is  station  azimuth.  S ,  P  and  B  are  source,  path  and  station 
terms  respectively.  S  and  B  depend  on  the  medium  structures  and  properties  at  the  source  and  receiver  locations, 
which  are  generally  different  in  a  heterogeneous  Earth.  The  path-dependent  term  P  has  the  form 

exp  -J  yR(coJ)dl 

- ,  (2) 

y  k(o))r0  sin  A 

where  A  is  the  epicentral  distance;  r0  is  Earth’s  radius;  U  is  the  group  velocity  of  the  surface  wave;  k  is  the 
wavenumber;  QR  is  the  quality  factor;  and  yR  is  the  attenuation  coefficient.  The  integral  is  taken  along  the  great- 
circle  path  between  the  source  and  the  receiver.  It  represents  the  average  attenuation  of  the  path. 

To  obtain  the  average  attenuation  from  measured  amplitudes,  we  calculated  theoretical  source  and  receiver  terms  S 
and  B  using  the  3D  velocity  model  CUB2  (Shapiro  and  Ritzwoller,  2002)  and  CMT  solutions.  We  then  removed  the 
source  and  receiver  terms  calculated  for  the  specific  source  and  receiver  locations  from  measured  amplitudes. 
Source-  and  receiver-corrected  amplitudes  were  used  to  estimate  the  average  attenuation. 


exp 


-co 


s 


dl 


2U((o,l)QR  {(0,1) 


Jk(co)r0  sin  A 
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18  s 


20  s 


Figure  4.  Path  coverage  of  single-station  amplitude  measurements  at  different  periods. 

Figure  5  plots  the  logarithm  of  corrected  amplitudes,  indicative  of  the  average  attenuation,  as  a  function  of  epicentral 
distance  at  12,  16,  and  20  seconds.  The  offset  from  zero  of  the  least-squares-fitted  line  at  zero  distance  indicates  that 
there  is  a  possible  bias  in  CMT  source  parameters,  consistent  with  the  observation  of  Yang  et  ah,  (2004).  In 
addition,  some  of  the  estimated  attenuation  coefficients  are  negative,  which  is  not  physical.  There  are  several 
possible  explanations  for  obtaining  negative  attenuation  coefficients.  They  include  un-modeled  site  response, 
remaining  contamination  in  the  measured  amplitudes,  errors  in  the  source  parameter  and  inadequate  description  of 
the  wave  propagation  by  the  ray  theory  (Equation  2).  We  then  removed  the  amplitudes  that  resulted  in  negative 
average  attenuation  coefficients  from  the  subsequent  tomographic  inversion. 

Attenuation-Coefficient  Tomographic  Inversion 

We  used  the  selected  source-  and  receiver-corrected  amplitudes  in  the  attenuation-coefficient  tomographic  inversion 
to  obtain  attenuation-coefficient  models  for  Asia  and  surrounding  regions  at  periods  between  12  and  22  seconds.  We 
used  a  modified  version  of  the  inversion  algorithm  described  by  Barmin  et  al.  (2001).  The  algorithm  inverts  for  the 
attenuation-coefficient  models  by  minimizing  the  functional 

^)  =  ^Z  (<4-)2>  (3) 

ij 

at  period  T,  where  N  is  the  total  number  of  paths  and  dy  is  the  amplitude  residual  between  source  i  and  station  j: 

d,j  =  qfs  -q°g  =  J  m(r)dl  +  In (SM,)  (4) 

V 

where  qfs  is  the  value  of  the  integral  in  Equation  2  from  measured  amplitude,  is  the  value  of  the  integral  from 
the  attenuation  model  that  we  invert  for,  and  m(r)  is  the  model  perturbation  from  a  reference  model  m0( r).  The 
additional  unknown  parameter  8Mt  is  included  to  account  for  the  possible  CMT  moment  bias. 
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12s 


Epicentral  Distance,  km 

Figure  5.  Source-  and  receiver-corrected  single-station  spectral  amplitudes  plotted  against  epicentral 

distance.  The  red  line  is  the  least  squares  fit  of  the  data.  Green  dots  are  data  rejected  from  the 
tomographic  inversion. 

We  minimized  F(T)  using  least  squares  and  several  damping  parameters  described  in  Barmin  et  al.  (2001)  with  an 
additional  damping  parameter  for  SMt.  Equations  3  and  4  were  used  in  the  inversion  of  18-s  amplitudes  to  obtain 
SMt.  The  resulting  M0i  =  M0  +  SMt ,  where  M0  is  the  CMT  scalar  moment,  was  then  used  to  correct  the  amplitudes  for 
the  source  term  at  other  periods,  and  SMt  was  removed  from  the  inversions.  Numerous  experiments  with  different 
values  of  damping  parameters  yielded  the  optimal  inversion  results.  We  also  experimented  with  inversions  in  which 
we  included  for  all  the  periods  (“free”  scalar-moment  inversion.)  The  results  are  similar.  The  resulting 
tomographic  models  of  attenuation  coefficients  jR  for  Asia  and  surrounding  regions  are  shown  in  Figure  6.  The 
attenuation  models  are  consistent  with  geological  features  of  Asia.  We  observe  low  attenuations  in  stable  regions 
such  as  eastern  Europe,  the  Siberian  platforms,  the  Indian  Shield,  the  Arabian  platform,  the  Yangtze  Craton,  and 
others.  High  attenuation  is  observed  in  tectonically  active  regions  such  as  the  Himalayas,  the  Tian  Shan,  Pamir,  and 
Zagros  mountains. 

We  estimated  variance  reductions  achieved  with  our  tomographic  models  by  comparing  their  residual  statistics  with 
those  of  homogeneous  models.  For  periods  between  12  and  22  seconds,  variance  reduction  is  between  30%  and  40% 
(Figure  7).  Including  a  moment  correction  term  SMt  in  the  inversion  did  further  reduce  the  variance. 

Calibration  of  Russell’s  (2006)  Ms  Formula  for  Asia 

Russell  (2006)  describes  a  new  formulation  of  surface-wave  magnitude  using  time-domain  and  Butterworth 
band-pass  filtered  amplitude.  The  new  Ms  is  applicable  for  amplitudes  measured  at  arbitrary  periods  and  is  not 
affected  by  the  dispersion  characteristics,  including  Airy-phase  anomalies,  of  the  surface  waves  because  of  the  band¬ 
pass  filtering.  The  ability  of  the  Russell  (2006)  Ms  to  use  amplitudes  at  variable  periods  is  important  in  calculating 
regional  surface-wave  magnitude  because  at  regional  distances  in  continental  regions,  surface  waves  usually  peak  at 
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Figure  6.  Tomographic  models  of  attenuation  coefficients  across  Asia  and  surrounding  regions.  Grey  color 

corresponds  to  areas  where  the  path  density  is  less  than  20  paths  across  an  equatorial  cell  of  2°  x  2°. 
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Figure  7.  Variance  reduction  of  tomographic  attenuation-coefficient  models  compared  with  homogeneous 

models.  In  “fixed-scalar-moment”  inversions,  SMi  obtained  from  the  inversion  of  18-sec  amplitudes 
was  used  to  correct  the  amplitudes  for  the  source  terms  at  other  periods.  In  “free-scalar-moment” 
inversions,  SM{  was  included  in  inversions  of  amplitudes  at  all  periods. 

periods  that  are  shorter  than  20  seconds,  at  which  traditional  surface-wave  magnitudes  are  defined.  The  Russell 
(2006)  Ms  is  defined  as 

Msm  =  log10  (ab)  +  ^  log10  (sin  A)  +  Batt  A  -  0.66  logJ-y] -  log10  {fc)  +  Cb  (5) 


and 


rV a 


where  ab  is  the  time-domain,  Butterworth  band-pass  filtered  amplitude;  A  is  the  epicentral  distance  in  degrees;  Batt  is 
related  to  the  attenuation  of  the  surface  wave;  T  is  the  period  of  the  amplitude;  and  fc  is  the  half  width  of  the 
Butterworth- filter  pass  band.  Gmin  is  the  minimum  value  of  G  where 


G  - 


U 


nb„ 


dU 


dT 


K 


(6) 


U  and  dU/dT  are  the  surface-wave  group  velocity  and  its  derivative  with  respect  to  a  period  for  the  region  of 
interest;  bn  is  related  to  the  order  of  the  Butterworth  filter;  and  at  is  the  degree-to-km  converting  factor.  For  typical 
continental  regions,  Russell  (2006)  found  that  Gmm  =  0.6  is  adequate.  Cb  in  Equation  5  also  depends  on  G. 

To  make  the  new  Ms  consistent  with  traditional  formulae  such  as  Rezapour  and  Pearce’s  (1998)  Ms ,  Russell  (2006) 
equated  Equation  5  with  traditional  formulae  and  derived  formulae  for  Batt  and  Cb: 


Batt  ~ 


Tov'8 


\TJ 


and  Cb-C  +  logj 


(Go. 

V400 


(7) 


In  the  above  formulae,  C  is  the  constant  in  traditional  Ms  formulae.  G0  is  G  at  20  seconds.  Adopting  the  formulation 
of  von  Seggem  (1977),  Russell  (2006)  obtained  a  B0  value  of  0.003 1  and  a  C  value  of  2.2  for  ab  in  nanometers. 
Using  nominal  values  of  U  and  dU/dT  at  20  seconds  for  continental  regions,  Russell  (2006)  derived  a  Cb  value  of 
-0.42  for  a  6th-order  Butterworth  filter. 

Bonner  et  al.  (2006)  developed  a  new  amplitude  measuring  technique  to  be  used  along  with  the  Russell  (2006)  Ms 
formula.  They  first  band-pass  filter  the  seismogram  between  8  and  25  seconds  using  a  bank  of  Butterworth  filters. 
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From  among  the  filtered  data,  they  then  find  the  time-domain  amplitude  that  maximizes  the  ratio  ab/fc,  thus  giving 
the  maximum  Ms  (Equation  5).  The  technique  yields  station  magnitudes  with  significantly  reduced  scatter  (Bonner 
et  ah,  2006). 

To  calibrate  Russell  (2006)  Ms  for  Asia,  we  collected  seismograms  from  98  events  that  occurred  from  2002  to  2003 
in  Asia.  We  used  Bonner  et  al.  ’s  (2006)  technique  to  measure  the  surface-wave  amplitudes.  We  then  calculated 
Russell  (2006)  Ms  with  Equation  5  using  the  nominal  parameter  values  derived  by  Russell  (2006).  We  also 
calculated  Prague  Ms  (Vanek  et  al.,  1962)  and  Rezapour  and  Pearce  (1998)  Ms  using  amplitudes  with  periods  in  the 
band  from  18  to  22  seconds.  The  Rezapour  and  Pearce  (1998)  Ms  is  currently  used  at  the  International  Data  Centre 
(IDC)  as  the  standard  Ms  measurement.  The  Russell  (2006)  Ms  that  we  calculated  using  the  nominal  parameter 
values  with  this  dataset  is,  on  average,  about  0.12  magnitude  unit  (mu)  larger  than  the  Rezapour  and  Pearce  (1998) 
Ms  calculated  with  the  same  dataset.  The  mean  station-magnitude  variance  of  the  Russell  (2006)  Ms  is,  on  the  other 
hand,  19%  smaller  than  that  of  the  Rezapour  and  Pearce  (1998)  Ms  and  58%  smaller  than  that  of  the  Prague  Ms , 
consistent  with  the  conclusion  of  Bonner  et  al.  (2006). 

To  derive  the  parameter  Cb  in  Equation  7  that  is  specific  to  Asia,  we  utilized  the  surface-wave  group-velocity 
models  developed  by  Levshin  et  al.  (2003).  We  first  converted  velocity  to  slowness.  We  then  took  the  arithmetic 
mean  of  the  slowness  for  all  model  nodes  as  the  average  slowness  for  the  whole  region.  The  average  group  velocity 
and  its  derivative  (Equation  6)  were  then  derived  from  the  average  slowness.  From  the  average  group  velocity  and 
its  derivative  at  20  seconds,  we  obtained  a  G0  value  of  0.72.  Using  this  G0  in  Equation  7,  we  obtained  an  Asia- 
specific  Cb  of  -0.54.  This  value  is  0.12  mu  smaller  than  the  nominal  value  of  -0.42  derived  by  Russell  (2006).  If  we 
use  this  value  in  Equation  5  to  calculate  the  Russell  (2006)  Ms  with  the  dataset  that  we  collected,  the  average  Ms 
will  be  the  same  as  the  average  of  Rezapour  and  Pearce  (1998)  Ms  from  the  same  dataset. 

Another  parameter  in  the  Russell  (2006)  Ms  formula,  which  is  also  region-dependent,  is  Gmin  in  Equation  5.  We 
found  Gm in  for  a  particular  path  from  among  the  G  values  along  the  path  calculated  from  the  group-velocity  model  of 
Asia.  We  then  use  the  smaller  of  Gm in  and  0.6,  a  value  suggested  by  Russell  (2006),  to  define  the  bandwidth  fc  of  the 
Butterworth  filters  used  in  filtering  the  data  and  in  the  Ms  calculation.  The  Russell  (2006)  Ms  calculated  using  the 
Asia- specific  Gmin  has  an  even  smaller  station-magnitude  variance.  The  mean  station-magnitude  variance  of  the 
Russell  (2006)  Ms  is  now  24%  smaller  than  that  of  the  Rezapour  and  Pearce  (1998)  Ms.  The  Russell  (2006)  Ms  itself 
is  increased  by  0.04  mu  on  average. 

CONCLUSIONS  AND  RECOMMENDATIONS 


We  developed  a  surface-wave  amplitude-measuring  tool  to  make  reliable  amplitude  measurements.  The  tool 
incorporates  techniques  such  as  phase-match  filtering  and  backazimuth  calculation  in  order  to  reduce  the  noise 
contamination  of  fundamental  Rayleigh- wave  amplitudes.  With  the  tool,  we  made  surface-wave  spectral-amplitude 
measurements  for  Asia  and  surrounding  regions  in  the  frequency  band  between  12  and  22  seconds.  Using  the 
amplitude  measurements,  we  conducted  tomographic  inversions  and  developed  2D  surface-wave  attenuation- 
coefficient  models  from  12  to  22  seconds  for  the  region. 

We  calibrated  the  Russell  (2006)  surface-wave  magnitude  Ms  using  the  2D  group-velocity  model  for  Asia  (Levshin 
et  al.,  2003).  Russell  (2006)  Ms  calculated  with  Asia-specific  parameters  reduces  the  station-magnitude  variance. 
The  average  magnitude  is  more  consistent  with  the  Rezapour  and  Pearce  (1998)  Ms  used  by  IDC.  We  plan  to 
conduct  further  calibration  study  of  Russell  (2006)  Ms  using  the  attenuation  models  that  we  developed. 
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